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KENDALL RANDOM WALK, WILLIAMSON 
TRANSFORM AND THE CORRESPONDING 
WIENER-HOPF FACTORIZATION 

B.H. JASIULIS-GOLDYN 1 AND J.K. MISIEWICZ 2 


Abstract. The paper gives some properties of hitting times and 
an analogue of the Wiener-Hopf factorization for the Kendall ran¬ 
dom walk. We show also that the Williamson transform is the best 
tool for problems connected with the Kendall generalized convolu¬ 
tion. 


1. Introduction 

In this paper we study positive and negative excursions for the Kendall 
random walk {X n : n £ No} which can be defined by the following 
recurrence construction: 

Definition 1.1. The stochastic process {X n : n £ No} is a discrete 
time Kendall random walk with parameter a > 0 and step distribution 
v if there exist 

1. (Y k ) i.i.d. random variables with distribution u, 

2. (6c) i-i-d. random variables with uniform distribution on [0,1], 

3. (9k) i.i.d. random variables with symmetric Pareto distribution with 
density TT 2a (dy) = a|t/| _2 “ _1 l[i i00 )(|?/|) dy, 

4. sequences (Yk), (6c) and (6k) are independent 
such that 


x 0 1 , X\ Tf, X n+ 1 Mi n j r \r n j r \ [I(6i ‘n Qn+ 1) T ^ri+ ll(6i dn+l)] i 
where 9 n+ \ and M n+ \ are independent, 

M n+ 1 = max{|A n |, |K„+i|}, m n+1 = min{|A n |, |K n+ i|}, g n +i = 

M n +1 


and 

r n +i = {sgn(u) : max{|A n |, |T n+ i|} = |tt|}. 
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The Kendall random walk is an example of discrete time Levy ran¬ 
dom walk under generalized convolution studied in [Ij. Some basic 
properties of this particular process are described in [3], 

In the second section we present the definition and main properties of 
the Williamson transform together with unexpectedly simple inverse 
formula. We give also the definition of the Kendall convolution ([3]) 
which is an example of generalized convolution in the sense defined by 
Urbanik M. This section is written only for the completeness of the 
paper, since all presented facts are known, but rather difficult to fold 
in the literature. 

In the third section we recall another definition of Kendall random walk 
following the multidimensional distributions approach given in HI. By 
this construction it is clear that the Kendall random walk is a Markov 
process. Then we study hitting times, positive and negative excursions 
for this process. The main result of the paper is the analog of the 
Wiener-Hopf factorization. 

For simplicity we use notation T a for the rescaling operator (dilatation) 
defined by (T a X)(A) = X(A/a ) for every Borel set A when a ^ 0, and 
T 0 X = ho- Li the whole paper parameter a > 0 is fixed. 

2. Williamson transform and Kendall convolution 

Let us start from reminding the Williamson transform (see ra and its 
basic properties. In this paper we apply it to ^-symmetric probability 
measures on R, but equivalently it can be considered on the set "Re¬ 
measures on [0, oo). 

Definition 2.1. By Williamson transform we understand operation 
is —> v given by 

v(t) — / (1 — \xt\ a ) + is(dx), iseV s , 

Jr 

where a + = a if a ^ 0 and a + = 0 otherwise. 

For the convenience we use the following notation: 

*(*) = ( 1-1*1“)+, G(t) = v{l/t). 

The next lemma is almost evident and well known. It contains the 
inverse of Williamson transform which is surprisingly simple. 

Lemma 2.1. The correspondence between a measure v e R and its 
Williamson transform is 1 — 1. Moreover, denoting by F the cumulative 
distribution function of is, i/({0}) = 0, we have 

p(+) — f TX i a (G(t) + 1) + tG'(t)] if t > 0; 

* K ) ~ \ 1 - F(-t) if t < 0. 

except for the countable many lei. 
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Proof. Of course the assumption z/({0}) = 0 is only technical simpli¬ 
fication. Notice first, that for t > 0 integrating by parts we obtain 


G{t) = (l-\x/t\ a ),dF(x) 


= 2 


(l — x a t a )F(x)\ t Q + at a x a 1 F(x)dx 


= — 1 + 2 at “ / x a F(x)dx. 


Since 


and 


lim — 

h— >-o+ h 


lim — 

h-> o- h 


ft+h 


x a F(x)dx — x a F(x)dx 


IJ o 


ft+h 


x a F(x)dx — / x a F{x)dx 


= t a ~ l F(t+), 


= t a ~ l F(t-) 


then, except for the points of jumps for F we can differentiate the 
equality 

2a [ x a - 1 F(x)dx = t a (G(t) + 1) 


and obtain 


2 at a ~ l F{t) = at 0 - 1 ( G(t ) + 1) + t a G'(t). * 

In the following we use the notation: S x = ^S x + \d- x and Tt 2 a ( dy ) = 
a\y\ — 2q! ~ 1 i[ i oo)(I?/1) dy is the symmetric Pareto distribution. 


Definition 2.2. By the Kendall convolution we understand operation 
A a : Vg —* V s defined for discrete measures by 

S x A a S y := T m (g a n2a + (1 - Q a )S 1) 

where M = max{|x|, |r/|}, m = min{|x|, \y\}, Q = m/M. The extension 
of A a to the whole V s is given by 


z'i A a u 2 (A) 


8 X A a 8 y (A) ui(dx)u 2 (dy). 


It is easy to see that the operation A a is symmetric, associative, com¬ 
mutative and 

• v A a = v for each v G V s \ 

• (pv 1 + (1 - p)v 2) A a O = p (zq A a v) + (1 —p) (z / 2 A a v) for each 
p G [0, 1 ] and each u, u\ , v 2 G V s . 

• if \ n -A A and v n -A v then (A n A a u n ) -A (A A a u), where —> 
denotes the weak convergence; 

• T a (y x A a 1+2) = {T a u 1 ) A a (T a v 2 ) for each v u v 2 G V s . 
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The next proposition, which is only a reformulation of a known prop¬ 
erty of Williamson transform, shows that this transform is playing the 
same role for Kendall convolution as the Fourier transform for classical 
convolution. 


Proposition 2.2. Let ui,u 2 G V s be probability measures with the 
Williamson transforms di,u 2 . Then 

/ '&{xt)(y i A a is 2 )(dx) = v x (t)v 2 [t). 

Jr 

Proof. Notice first that for M = max{|x|, |y|}, m = min{|x|, |y|},: 
g — m/M we have 

f ty(ut)(5 x A q 5y)(du) = ( ty(Mut) (<5 e A a <5i) (du) 


= g a 2 / 'ff(Mut) 


a du 


u 


\2a+l 




= g a (1 - \Mt\ a f l[_i,ij(Mt) + (1 - g a ) (1 - | Mt\ a ) 1 [— 1 , 1 ] (Mt) 

= (1 - \Mt\ a ) (1 - \mt\ a ) l[_i,i ](Mt) = ^(xt)^(yt) 

Now, by the definition of Kendall convolution, we see that 

/ \&(xt)(i/i A a u 2 )(dx) — / ^(ut) (S x A a S y ) (du) vi(dx)v 2 (dy) 

Jr Jr 3 

= f ^(xt)^(yt)u 1 (dx)u 2 (dy) = (t)u 2 (t), 

Jr 2 

which was to be shown. 4k 


Proposition 2.3. There exists a sequence of positive numbers (c n ) 
such that T Cri 5f aU converges weakly to a measure cr ^ 5 0 (here u Aan = 
v A a ■ ■ ■ A a v denotes the Kendall convolution of n measures u). 

Proof. Since the Williamson transform of S\ is equal \l/(f) then putting 
c n = rC 1 ! 01 we obtain 


V(xt)T c J? an (dx) = / ^(c n xt)8f an (dx) 


= (1 - \c n t\ a ) 


a\n 
+ 


Tl“ if 


n —>• oo. 


The function e~^ a is symmetric and monotonically decreasing as a 
function of t a on the positive half line, thus there exists a measure 
a G V s such that 


e-W 



xt\ a ) + cr(dx). 


According to Lemma 2.1 we can calculate the cumulative distribution 
function F of the measure a obtaining F(0) — J and 


F(t) 


|(1 Ft~ a + e t °‘)e~ t “ if t> 0; 
1 — F(—t ) if t < 0. 


4k 
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Remark. Notice that properties of Kendall convolution together with 
the one described in the last proposition show that the Kendall con¬ 
volution is an example of generalized convolutions in the sense defined 
by K. Urbanik. More about generalized convolutions one can find in 

EM 12 El ME]. 

As a simple consequence of Lemma 2.1 and Proposition 2.2 we obtain 
the following fact: 


Proposition 2.4. Let v G V. For each natural number n ^ 2 the 
cumulative distribution function F n of the measure i/ Aa11 is equal 


F n (t) = — [a(G(t) n + 1) + tnGity^G'it)] , t > 0, 

and F{t ) = 1 — F(—t ) for t < 0, where G(t) = v(l/t). 

Examples. 

1. Let v be uniform distribution on [—1,1], Then 


G(t) = 




OL - 1 “ 1 
1 1 


OL 1 


l(l,oo)(|t|)> 


a(i) = 2 + 2 sgn(t) (^ T )“ (t + 1 ) r-'ivMt) 


+ \ Sgn(f) (l - 


t~ a \«-i 


ol T 1 


n t 
1 + - + 


a a + 1 V t 


n 


l) ^l[l,oo)(^). 


2. Let u — 8i \8\ + |h_i, thus G(t) = (1 — \t\ “)+. For each 
natural number n ^ 2 the probability measure has the density 

dF n (t) = a ’* ( t "~ 1 11 (1 - |t|-“)— 2 (|(|) dt 


3. Let v = p5\ + (1 — p)tt p , where p G (0,1] and tt p is the symmetric 
Pareto distribution with density TT p (dy) = a|y| _p_1 l[i,oo)(|z/|) dy. For 
each natural number n > 2 we have: 


F n(t) = ^ + isgn(t) 


l-°( 1 -P)| tr + P( 1 -°)| t |-c 1 "- 1 


x 


1 + 


a — p a — p 

(1 - p)(np - a) p( 1 - a)(n - 1) 


a — p 


-i*r 


a — p 


\t\ “ l[l,oo)(M). 


4. Let v = 7 r 2a for a G (0,1]. Since hi A a hi = 7r 2Q then using 
Example 2 we arrive at: 


dF n (t ) = 


an(2n — 1) , , ,_ a \ 2 (n-i). 


\t\ 


2a+l 
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3. Kendall random walk 

The direct construction of the symmetric Kendall random walk {X n : 
n G N 0 } based on the sequence (Y k ) of i.i.d. unit steps with distribu¬ 
tion v we have already presented in Definition 1.1. Slightly different 
direct construction was given in [2]. We recall here one more definition 
of the Kendall random walk following [I], where only multidimensional 
distributions of this process are considered. This approach is more con¬ 
venient in studying positive and negative excursions, which we consider 
in this section. 

Definition 3.1. The Kendall random walk or a discrete time Levy 
process under the Kendall generalized distribution is the Markov process 
{X n : n G N 0 } with X Q = 0 and the transition probabilities 

P n (x, A) = P {X n+k G A\X k = x} = 6 X A a z/ A “ n , n, k G N, 

where measure v G? s is called the step distribution. 

For any real a we introduce the first hitting times of the half lines 
(a, oo) and (—oo, a) for random walk {X n : n G N 0 }: 

T+ = min{n ^ 1 : X n > a}, r~ = min{n ^ 1 : X n < a} 

with the convention min0 = oo. We want to find the joint distribution 
of the random vector (tq 1 ", X t +). In order to attain this we need some 
lemmas: 


Lemma 3.1. 


1 / xy 

5y(0,t) 2 \ 


1 r 

2 




X 




y 




x 




y 


L ii a: i<ti’/i<fi) 


5 x A a u(0,t) = Pi{x,[0,t)) 








Proof. With the notation a = min{|a:|, |y|}, b = max{|x|, ||/|} and 
g = (a/b) a we have 


S x A a 5 y (0, t) = T b ((1 - g)5i + qtt 2q ) (0, t ) 


= ^(1 - £010 < t) + ^gT b n 2a (0,t) 

= - g)l(b < t) + ^Ql(b < t) I 


pt/6 2a 


1 
2 

= 0" 


0 - 0 ^ 0 *->^ 0 -©><<> 
l(|x| < t, \y\ < t). 


s 2a+l 

fab\ a 


ds 


xy 

t 2 
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The final form trivially follows from the identity 1 — ab = 1 — a + 1 
b — (1 — a) (1 — b). The second formula follows from 


AQ ^(0,t) / Aq 6y (0, t ) h'^dy'j . 


X 


The next lemma is crucial for further considerations. 


Lemma 3.2. Let { X n : n G No} be the Kendall random walk. Then 


$n(t) ■= P{Xi < 0,... X n _i < 0, 0 < x n < t} 


= - G(tr~L 2n lF(t) - - - (n - l)G(t) 


Proof. For k 

k = 2 


= 1 we have $i (t) = P{0 < Ad < t} — F(t) — |. For 


/ 0 r*t r*0 

/ Pi(xi,dx 2 )Pi(0,dx 1 ) = Pi(x l ,(0,t))u(dx l ) 

-oo J 0 J — oo 

1 />0 , 


’ —t L 


2T(xi/t) ( F(t) - - ) + G{t) - T(xi /t)G(t) v(dx i) 


4 L 


4G(t)(F(t) 


G(tY 


In order to calculate $3 notice first that by Proposition 12.21 


V(y/t) (S x A a v) ( dy ) = 4>(x/t)G(t) 


and applying the second formula from Lemma 3.1 to the measure 
Pi(xi, •) = S X1 Aq v we have 


1*0 r*t pO 

/ / Pi(x 2 ,dx 3 )P 1 (x 1 ,dx 2 ) = / (5 ) (0, t) (5 x\ A a v) (dx 2 ) 

' —00 Jo J —00 

1 


1 —t 


2*(x 2 /t) \F(t) - - ) + G(t) - Or 2 /t)G(t) 


(< 5 Xl Aq v) (dx 2 ) 


1 

2 

1 

4 


^(xi/f)G(f)(F(f) - +G(t)S xl Aq i/( 0, t) - ^(xi/t)G(t) 2 


4V( Xl /t)G(t)(F(t) - - ) + G(t) 2 - 24>( Xl /t)G(tf 
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Consequently 

/*0 rO 

**(t) = I I 

l rO r 


Pi(x 2 , dx 3 )Pi(xi, dx 2 )Pi{0 , dx i) 


—oo J —oo J 0 


4 

1 

4 


'-t 


4^/(xi/t)G(t)(F(t) — ^ ) + G(t) 2 — 2'fr(x 1 /t)G(t) 2 \ u(dx i) 


2 G(tf(F(t) - i) + G(tf(F(t) - i) - G(tf 
6 G(t) 2 (F(t) - i) - 20(0= 


1 

23 


Simple, but laborious, application of mathematical induction ends the 
proof. X 

Lemma 3.3. The random variable (and, by symmetry of the Kendall 
random walk, also variable Tq ~) has geometric distribution P(tq = k) = 


2 k i 


k = 1 , 2 , • • • it has the followmg moments generating function: 


E.s T ° = 


1 - £’ 
- 1 - 2 


0 < s < 2 . 


Proof. Notice that lirn^oo G{t) = 1 since lim^^ ^(l/f) = 1. Now it 
is enough to apply Lemma 3.1: 

1 


P( T o =k) = P {X ] ^ 0 , • • • , X k _i ^ 0 , X k > 0 } = <h fc (oo) = 


2 k ' 


X 


Lemma 3.4. 


P \X + <t\ = 


4 Fit) - 2 - G(tf 


(2 - G(t)f ■ 
Proof, ft is enough to apply Lemma 3.1: 

OO OO 

p W < *} = £ P { < *> r 0 = k } = £ 


k =1 


k =1 


E ^ \ 2n ( F w - o) - («- 


n=l 


(2 - G(f))* 


4F(t) - 2 - G(f) s 


X 


Definition 3.2. Let {X n : n G N 0 } be the Kendall random walk in¬ 
dependent of the random variable N s / 2 with the geometric distribution 

P{N s/2 = k} = (f) fc_1 (l ** |), 0 ^ s ^ 1, for k = 1,2,- By the 

geometric-Kendall random variable we understand the following 


Z s /2 - E Xk l 


{K/2=k}- 


k =1 
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Notice that ET (X k /t) = G(t) k by Proposition 2.2, thus we have 


E *{Z s/2 /t) 




k =1 




gW(i-f) 

l-f G(t) ' 


We are ready now to prove the main result of the paper: the Wiener- 
Hopf factorization for the Kendall random walk. This gives the Laplace- 
Williamson transform Es T o T (uX T +), thus consequently, also the joint 
distribution of ('r 0 l ",A" r +). 


Theorem 3.5. Let {X n : n G N 0 } be a random walk under the Kendall 
convolution with unit step X\ rs_/ F such that F( 0) = Then 


E s T ° T 


K„ + ) 


I (1 - l)G(l/u) 
1 - f ' 1 - §G(l/w) 


Es T ° + ET ( uZ s/2 ) . 


Proof. Let H(s,u) := Es T o T (juX r +J. Then 

oo 

OO /*00 p / OO 

H(s,u ) = / X &(ux)d$ k (x) = I ^(ux) d\ ^ s k $ 

k =l •/ k fe= i 


fe(^) 


- 2 J - (* ~ !)G(a;) 


= / T(-ua;)d 


s(F(x) - |) - ±s 2 G(a;) 2 


= aw 


(l — |sG(x)) 
1/ “ z , Q _i g(E(a;) - |) - \s 2 G(x) 

(l — 


dx. 


In these calculations we write d(f(s,x)) x for the differential of the 
function f(s,x ) with respect to x. To get the last expression we need 
to know that G(0) = 0. Now we calculate G'(x). Since 


G(x) = / T (—\ v(dz) = 2 / (l — — \ u(dz), 

Jr Jo ' x / 

then integrating by parts we obtain 

pX 

(G(x) + 1 ) x a = 2a / z a ~ 1 F(z)dz. 

Jo 

Now it is enough to differentiate both sides of this equality to obtain 
G'(x) = 2ax~ 1 (f{x) - - - -G(x) 
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We see that 

H(s, u ) = u c 


j-l/u 

i 

TT 

s 

i_ 

L 

_1-| G{x)_ 


dx = 


\G(l/u) 


1 - §G( 1 /m)' 


The same method as in the proof of the last theorem we are using in 
the next lemma. 

Lemma 3.6. 

EX a + = 2E\Y 1 \ a 

T o 11 

Proof. In the following calculations we take R > 0: 

'F{x)-\-\G(x) 2 ' 


EX% = 


x a d 


= lim 

R —>oo 


= lim 

R—yoo 


= lim 




X 


F(x) 


(1 - \G{x)f 

| - \G{xf 


(1 - \G{x)) 2 
F(x) - | - ±G(x) R 


r r R 

- x a \G(x ) ■ 

/ 

dx 

o Jo 




(i - ±G(x)y 


x 


R->°° L (1 — kG(x)) 2 - 


u(du) — 


PX 

Jo 

l - 

K 

SL 

_1 


\xJ . 

J . 


R 


= lim 


R^oc [(1 - kG(x)) 2 J 0 


u a is(du ) 


R 


= 4 


u a is(du), 


J o 


which was to be shown. Notice that the assumption E|Yi|“ < oo is 
irrelevant in these calculations. X 
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